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Abstract 

We present an analytically tractable variant of the voter model that provides a quantitatively 
accurate description of /3-diversity (two-point correlation function) in two tropical forests. The 
model exhibits novel scaling behavior that leads to links between ecological measures such as 
relative species abundance and the species area relationship. 

PACS numbers: 87.23. Cc,87.23.-n,05.40.-a,05.90.+m 
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An ecological community represents a formidable many-body problem - one has an inter- 
acting many body system with imperfectly known interactions and a wide range of spatial 
and temporal scales. In tropical forests across the globe, ecologists recently have been able to 
measure certain quantities such as the distribution of relative species abundance (RSA), the 
species area relationship (SAR), and the probability that two individuals drawn randomly 
from forests a specified distance apart belong to the same species (also called /^-diversity) . 
In this letter, we take a first step towards the development of an analytically tractable 
model that, despite its simplicity, leads to a remarkably accurate quantitative description of 
/3-diversity in two different tropical forests. It also indicates the existence of novel scaling 
behavior, revealing previously unexpected relationships between /3-diversity, RSA, and SAR. 
The model we study is a version of the well-known voter model 111 which has been applied 

nnn 

to a variety of situations in physics and ecology [2, y, 1^ • 

Quite generally, one may study an ecosystem in d spatial dimensions with the most 
common case corresponding to d = 2. We will consider a hyper-cubic lattice in d dimensions 
with each site representing a single individual and where the lattice spacing, a, is such that 
a'^ is of the order of the average volume (area in d = 2) occupied by an individual. At each 
time step an individual chosen at random is killed and replaced, with probability 1 — u with 
an offspring of one of its nearest neighbors or, with probability u, with an individual of a 
species, not already present in the system. This last process is called speciation and the 
parameter u is called the speciation rate. The case with z/ = is special (the i/ — >■ case 
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is distinct from u = 0) and has been thoroughly studied (see, for example, 
on a finite lattice results in a stationary state with just a single species (monodominance). 
The case u = 1, in which a new species is generated every time step, is trivial. Our focus 
is on < < 1, which is ecologically relevant and has been studied before j2[|^, but is not 
well understood. 

Let Fg be the probability that two randomly drawn individuals separated from each 
other by x at time t are of the same species (for simplicity, the system is assumed to be 
translationally invariant). The master equation for a community of size occupying an 
area A = a'^N may be written as [lol | 



where {cj} is the basis vector set and x 7^ 0. When a; = Fq = 1. 
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a 


b 


z' 


1 


0.5 


1.05 


0.03 


2 


0.87 


1.2 


0.3 


3 


0.95 


1.21 


0.5 



TABLE I: Scaling exponents for d = 1, 2, 3 determined from the scaling collapse of the normalized 
RSA and SAR plots (Figs. 1 and 2). 



In the continuum limit, Eq.Q becomes 

dF{r, t) 



dt 



(2) 



where r = F{r,t) = a^'^Fi, 7^ = 2dv/a'^ and the time is measured in units of r = 
cr^/ {2dN). The continuum limit is obtained on choosing 00, cr — > and — in such 

a way that a'^N and 7 approach constant values and r — > 0. The first term on the right 
hand side represents dispersal or diffusion. The second term is a decay term arising from the 
effects of speciation jj], whose coefficient, 7^, could generally be a function of r. The last 
term is a consequence of the fact that for the discrete case, Eq(^, at x = one necessarily 
has the same species by definition. The constant c is fixed such that j^^ d'^r F{r) = 1 where 
C is the cube of side a (the average nearest neighbor plant distance) centered in the origin 

H. 

The stationary solution of Eq.((2I) is 



F(r) 



sr 2 K2-d ('yr). 
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(3) 

We 



where K is the modified Bessel function of the second kind [1 
have carried out extensive simulations on a square lattice and on a hexagonal lattice with 
periodic boundary conditions in both cases and have verified that the results are in excellent 
accord with the analytic solution and independent of the microscopic lattice used. 

When d = 1 the above expression takes a simple form: F{r) ~ exp(— 7r). It is easy 
to verify that the stationary solution of the discrete equation (^Q) with (i = 1 is also an 



exponential function Fg = e where \ = a 



In- 



is the correlation length. 



In the small v limit A ~ 7 



l-^y(2-u)_ 

^1 \fv and one may identify a natural scaling variable Av'^l'^ . 
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A second scaling variable is identified by noting that 
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■JL, x,y \ i=l / 

provides a characteristic scale for ra, where the indicator I^^gis a random variable which takes 
the value 1 with probability F^^g and with probability (1 — -Fk-j/), 5* and represent the 
number of species and the number of individuals of the i-th species respectively and (. . .) is 
an ensemble average, or a time average over a long time period. Furthermore, using Eq.(^ 
and its Fourier Transform, one finds that ^^-^^ ~ / -F(r)d'^r ~ z/~" with a = min(l,(i/2). 
This results in a second natural scaling variable nz/". 

Within the context of the same model, we turn to an investigation of other quantities 
of ecological interest, notably the relative species abundance (RSA) and the species-area 
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relationship (SAR) [13|, [IJ] in the limit of small speciation rate, u, and large volume (area 
for d = 2) A of the system. Here we will investigate the case when i/ is small but the system 
is still far from the onset of monodominance {Av'^^'^ > !)• 

First, let us postulate a scaling form for the total number of species S'(z/, A) within an 
area A: 

S{u,A) = A''S{Aiy''/^). (4) 

The exponent z' is equal to the traditional species-area relationship exponent, z, only if S 
is constant. From our numerical simulation we find that this is true when A <^ z/"'^/^, while 
when A ^ u~^^'^ S{u^ A) is linear in A. 

Let us introduce a scaling form of the normalized RSA distribution of species, 0(n, A) = 
/(n, z/. A)/ S{u, A), where /(n, z/. A) is the number of species with n individuals: 

I/, A) = n-^0(nl/^ Au'^^^), (5) 

where b is an as yet undetermined exponent. By definition, 

^0(n,z/,A) = l (6) 

n 

and 

J2nHn,'^,A)^A/S{u,A) (7) 

n 

(for simplicity, we set a = 1 so that A = N). 
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FIG. 1: Left column: plots of the normalized RSA for d = 1,2,3 with = 
0.001,0.003,0.01,0.03,0.1 (the d = 2 plot also shows the results for u = 0.0001,0.0003,0.3); 
L = 200. Right column: plots of the data collapse yielding a measure of the exponents a and 
b in Table 1. 



The normalization condition Eq.® becomes 

POO 



(8) 



Assuming that y) does not diverge as x tends to zero (which is justified a posteriori by 
numerical simulations), one finds that 6 > 1 if (f){x,y) approaches a constant value (modulo 
multiplicative logarithmic corrections) and 6 = 1 if 0(x, y) as x — ^ 0. In both cases 
(f){0,y) is independent of y (also in accord with the results of computer simulations). The 
condition on the average population per species, Eq.©, leads to 



(9) 



where S'{y) = y^ ^' /S{y). Detailed simulations followed by a scaling collapse indicate that 
in all dimensions, b < 2. It then follows that the lower integration limit in the above 




FIG. 2: Left column: plots of the SAR for d = 1, 2, 3 with v = 0.001, 0.003, 0.01, 0.03, 0.1. Right 
column: plots of the data collapse yielding a measure of the exponent z' in Table 1. 

integral can be safely put to zero (a non-zero lower limit of integration would merely result 
in corrections to scaling) leading to the scaling relation 

ai2-b) = ^{l-z') (10) 

and 

POO 

S'{Au'^^^)= dxx^-''(^{x,Aiy'^/^). (11) 
Jo 

The linear dependence of S'(z/, A) on A arises from becoming independent of A in the limit 
A ^ z/"'^/^ (see Eq.Q). This then leads to the scaling function S{y) ~ y^~^' in Eq.((Z|). This 
also follows from noting that 5" approaches a constant value for large A (Eq.Q). 

We have carried out extensive simulations with hypercubic lattices of various sizes in 
d = 1, 2, 3. A series of simulations with fixed size and varying speciation rate was used for 
the determination of the normalized RSA {L = 200 for d = 2 and L = 100 for d = 1,3, 
L being the side of the hypercube used). Another series of simulations, varying both the 
speciation rate and L was carried out to determine the SAR curves. 5(1^, A) is the mean 
number of species in a simulation with speciation rate u on a. hypercubic lattice of size 
A = L'^. In order to carry out the collapse, we used the automated procedure described in 
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FIG. 3: Beta Diversity data along with the best fits using Eq. (|12j) for plots in (a) Panama 
{R = 46m, = 68m, 7^^^ = 210km and cq = 0.012) and in (b) Ecuador-Peru (Yasuni) {R = 
86m, 7q"^ = 69m, 7^"^ = 23,500km and cq = 19). The Janzen-Connell effect pushes conspecific 
individuals further away from each other and thus the probability function Fq declines more steeply 
within the zone of its operation than at larger distances as in the data. Note that for length scales 
much bigger than R, in a regime in which Janzen-Connell effects do not play a significant role, the 
beta diversity curve decays much slower with length in Yasuni than in Panama by nearly a factor 
of 100. This difference may be attributed to a very strong climatic gradient (in annual rainfall) 
across the isthmus in Panama, which causes rapid spatial changes in forest composition. In Yasuni, 
on the other hand, there are large areas with very similar climate. The experimental protocol and 
a description of the symbols are in 



1(t| . Applying this procedure to the data on the normahzed RSA obtained in our simulations, 
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we were able to obtain the values of the exponents a and h (see Table HI and FigH]). 

As expected, our extensive computer simulations indicate that (f) is only weakly dependent 
on Au'^/'^. In all dimensions the scaling exponents (see Table HJ) approximately satisfy Eq. ljlUp 
and h is in the interval (1,2). Figure |21 shows a collapse plot which confirms the scaling 
postulates above. The biggest deviation from our theory is found in the value of a in 
(i = 2, the upper critical dimension for diffusive processes and are suggestive of logarithmic 
corrections. Interestingly, our scaling relation seems to hold even for this case. 

We conclude with a generalization of our model which provides an excellent quantitative 
fit to the /3-diversity data of two tropical forests. The key idea is that the factor 7 in 
Eq.Q can, in principle, take on two distinct values (7 = 70, r < _R; 7 = 71, r ^ R, 
where i? is a characteristic length scale separating the two distinct regimes). Such an effect 
arises physically from now well-established density-dependent processes first hypothesized by 
Janzen jl^ and Connell in which the survival probability for an offspring is decreased 
in the vicinity of adults of the same species. Janzen and Connell postulated that this 
increased mortality rate of seeds and seedlings near adults arise from the presence of pests 
that are host-specific, i.e, specialized to that type of tree, and experimental evidence supports 
this conclusion The microscopic model that we consider is slightly modified from the 
previous version in that there is a non-zero probability that a new-born is immediately killed 
with a probability proportional to the number of conspecific adults within a circle of radius 
R centered at the site. 

The solution of Eq.(j2I) in d = 2 with two distinct values of 7 is 



CoKoijor) + Ci/o(7or), r < R 

(12) 

C2Ko{'yir), r > R 



Fn(r) 



where Iq and Kq are modified Bessel functions [12| and the constants ci and C2 are fixed by 
imposing the continuity of Fo{r) and its derivative at r = i? (-^o('") diverges as r ^ 00 and 
is therefore excluded in the region r > R). Figure |S1 shows that our theory leads to good 
fits of the data on /^-diversity for tree communities in both central Panama (top panel) and 
Ecuador- Peru [isj l. 
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